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1.0 INTRODUCTION 

The success of imaging spectrometry in mineralogic mapping of natural terrains 
(e.g., Clark et al., 1992) indicates that the technology can also be used to assess the 
environmental impact of human activities in certain instances. Specifically, this paper 
describes an investigation into the use of data from the Airborne Visible/Infrared Imaging 
Spectrometer (AVIRIS) for mapping the spread of, and assessing changes in, the 
mineralogic character of tailings from a major silver and base metal mining district. 

The area under investigation is the Coeur d'Alene river valley in northern Idaho. 
Mining has been going on in and around the towns of Kellogg and Wallace, Idaho since 
the 1880‘s. In the Kellogg-Smelterville Flats area, west of Kellogg, mine tailings were 
piled alongside the South Fork of the Coeur d*Alene River. Until the construction of 
tailings ponds in 1968 much of these waste materials were washed directly into the South 
Fork. The Kellogg-Smelterville area was declared an EPA Superfund site in 1983 and 
remediation efforts are currently underway. Recent studies (e.g., Horowitz et al., 1992) 
have demonstrated that sediments in the Coeur d’Alene River and in the northern part of 
Lake Coeur d’Alene, into which the river flows, are highly enriched in Ag, Cu, Fby Zn, 
Cd, Hg, As, and Sb. These trace metals have become aggregated in iron oxide and 
oxyhydroxide minerals and/or mineraloids. Reflectance spectra of iron-rich tailing 
materials are shown in Figure 1. Also shown in Figure 1 are spectra of hematite and 
goethite. The broad bandwidth and long band center (near 1 pm) of the Fe^"^ crystal-field 
band of the iron-rich sediment samples combined with the lack of features on the Fe^^-0^‘ 
charge transfer absoiption edge indicates that the ferric oxide and/or oxyhydroxide in these 
sediments is poorly crystalline to amorphous in character. Similar features are seen in 
poorly crystalline basaltic weathering products (e.g., palagonites). 

The problem of mapping and analyzing the downriver occurences of iron rich 
tailings in the Coeur d'Alene (CDA) River Valley using remotely sensed data is 
complicated by the full vegetation cover present in the area. Because exposures of rock 
and soil were sparse, the data processing techniques used in this study were sensitive to 
detecting materials at subpixel scales. The methods used included spectral mixture 
analysis (Adams et al., 1993) and a constrained energy minimization technique (Farrand 
and Harsanyi, 1994a). 

2.0 DATA 

The CDA River Valley is currently under investigation by the Spokane field 
office of the USGS (C. Smith, pers. comm.). In a cooperative endeavor between the 
USGS, the EPA and the Naval Research Laboratory’s HYDICE (Hyperspectral Digital 
Imagery Collection Experiment) Project Office an AVIRIS flight was scheduled over the 
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CDA River Valley. Data was collected on May 24, 1993 at 21:40 UTC. Incomplete 
coverage of the target area led to a reflight on September 17 at 19:40 UTC. An earlier 
paper (Farrand and Harsanyi, 1994b) presented a preliminary examination of two scenes 
from the May collection. The current paper considers three scenes from the September 
collection that cover an area extending down river from approximately Smelterville to 
Killamey Lake (see Figure 2). 

The reflectance of samples collected in the field were measured in the laboratory. 
A subset of the samples were measured at Brown University’s RELAB. The majority of 
the samples were measured on the USGS Reston Beckman (UV5240). 

The AVIRIS data were analyzed in a radiometrically calibrated fmnat and also in 
a format with DN equal to apparent surface reflectance. The conversion from radiance to 
reflectance was done using both surface reflectance data and the ATREM (Gao et al., 

1993) radiative transfer based approach. First the data were converted to parent surface 
reflectance through the use of ATREM which effectively removed atmospheric variability. 
This was followed by the application of a modified flat field (MFF) method (Farrand, 
1992) applied to the ATREM corrected data which acted to remove instrumental noise that 
is retained in the ATREM correction. 

3.0 ANALYSIS 

The AVIRIS radiance data were initially examined with spectral mixture 
analysis. An iterative approach utilizing several runs and successive examination of the 
resulting RMS error images (Adams et al., 1993) suggested four endmembers: shade, 
vegetation, an agricultural soil, and an iron oxide or oxyhydroxide rich soil. As was 
discussed above, the iron-rich soil is the most likely candidate for bearing heavy metals 
leached from the mine tailings. In the fraction image for this endmember, other materials 
besides the ferruginous soil (for instance, the pavement on 1-90) appeared bright on the 
fraction image. 

In order to cut down on the false alarm rate, a recently developed detectiem 
algorithm that has already shown good success in detecting distributed subpixel target 
materials, the constrained energy minimization (CEM) technique (Farrand and Harsanyi, 
1994b), was applied to the data. The key to the CEM technique is to determine a vector 
operator, w composed of weights (wi...w^ that suppresses the unknown and undesired 
background spectra while enhancing the target spectrum d. The summed pixel energy can 
be represented by a scalar value y^ . The CEM operator is defined by two constraints. 

The first constraint, for any given pixel, is to minimize the energy summed across the 
wavelength range (e.g., minimize y.). The second constraint is that when applied to the 
target spectrum, y, = 1, e.g., w^d=l. 

CEM was applied to subsections of the three AVIRIS scenes. The target 
spectrum for input to the CEM routine was derived from the data itself. Three pixels 
most like the BH-3 sample spectrum shown in Figure 1 were identified. The radiance 
spectra of these pixels were then extracted from a radiance data set that was subsectioned 
to 188 bands to exclude several of the shortest and longest wavelength bands, 
spectrometer overlap regions and atmospheric water absorptions at 1.4 and 1.9 pm. 

These radiance spectra were averaged and the average was used as input to the CEM 
routine. As might be expected, the resulting CEM component images showed that the 
iron-rich sediments are concentrated primarily along the CDA River. There were also 
known concentrations of tailings in a dump west of Smelterville. Weaker responses are 
returned from some fields in the CDA River Valley. 

The CEM routine identified 159 pixels in the three scenes that were most like 
the iron-rich sediments at Cataldo Flats. These pixels were then extracted from the 
ATREM+MFF corrected data cube and transformed via a principal components analysis 
(PCA). The PCA was applied in two parts. 88 VNIR channels were analyzed separately 
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from 45 SWIR channels. This was done because iron bearing minerals primarily have 
features in the VNIR and clay minerals primarily have features in the SWIR. By 
separating the analyses, the resulting intrinsic dimensionality of the data set could be kept 
to a manageable number of endmembers (e.g., 3 to 4). For purposes of brevity only the 
VNIR analysis will be discussed here. 

The first three principal components of the VNIR analysis were examined 
interactively using the "XGobi" data visualization program (Swayne et al., 1992). By 
rotating the data cloud, spectrally extreme endmember pixels could be identified. Two 
distinct linear trends revealed themselves outlining two ridges of a tetrahedron. Pixel 
spectra from the ends of these trends (i.e., the vertices of the tetrahedron) are shown in 
Figure 3. Although the PCA was done on only the first 88 channels, the full spectral 
range is shown in Figure 3 since, as it turns out, there is extra information to be gained 
from the SWIR region as well. Figure 3 reveals two similar af^iearing, but nonetheless 
distinct, sets of spectra. Both trends have one endmember with a strong "1 pm" Fe 
crystal field absorption feature and one with only a weak absorption at that wavelength. 
The spectra, t2eml and t2em2 appear to be wetter with stronger 1.4 and 1.9 pm 
abswption feauues. These spectra also display absorptions indicative of gypsum at 1.48, 
1.53 and longwards of 1.7 pm. The presence of the sulfate mineral, gypsum, lends 
further credence to the interpretation that the ferruginous sediments are contaminated by 
mine wastes. 

4.0 CONCLUSIONS 

These results provide encouragement that imaging spectrometer data can be used 
effectively both to track the spread of fluvially distributed mine tailings and to provide 
insight into changes in their mineralogic character. An enabling algorithm that aids 
significantly in this process is the constrained energy minimization technique which, in 
this instance, provided a unique identification of the ferruginous sediments associated with 
mine waste contamination. 
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figure 1. reflectance spectra of iron-rich sediments collected near the Cocur d'Alene River 

cewnpared against library spectra of goethite and hematite. 
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Figure 2. Sketch map of Coeur d'Alene River Valley in nonhem Idaho. The AVIRIS flightlincs arc 
outlined by the box. 
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Rgurc 3. Iron-rich sediment endmembers determined firom PCA of VNIR channels and 
interactive data analysis. 
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